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Abstract 


Numerical  solution  of  the  quasi-one-dimen- 
slonal  unsteady  compressible  flow  equations  applied 
to  flow  through  the  control  valve  and  slotted  air 
supply  duct  of  a  circulation  control  rotor  is 
described.  An  implicit  finite-difference  method 
is  used  for  the  numerical  solution.  Solutions  are 
compared  with  basic  experimental  results  obtained 
for  a  mock-up  of  a  circulation  control  rotor  and 
its  pneumatic  valving  system.  Details  of  the  ex¬ 
periment  are  presented  and  the  experimental  and 
numerical  results  are  discussed. 

NCHENCLAIURE 

A  Cross-sectional  area  of  duct 

Ae  Effective  expansion  area  at  valve  exit 

A^  Valve  area 

a  Jacobian  matrix  3F/3U 

b  Jacobian  matrix  3G/3U 

Cf  Duct  friction  factor 

C  Slot  discharge  coefficient 

Cr  Tip  discharge  coefficient 

tip 

cv  Specific  heat  at  constant  volume 

Dg  Hydraulic  diameter  of  duct 

e  Internal  energy  per  unit  volume 

F  Flux  time  a  area  vector 

G  Vector  of  nonhomogeneous  terms  in  flow 

equations 

h  Heat  transfer  coefficient 

Hp  Average  duct  height 

hp  Total  enthalpy  in  pressure  supply  plenum 

pi  minus  rotational  kinetic  energy  at  duct 
entrance 

k  Thermal  conductivity 

m  Mass  flow  rate  per  unit  area  in  duct 

«j  Mass  flow  rate  per  unit  area  computed 

from  isentroplc  flow  theory 
ma  Mass  flow  rate  per  unit  area  through  slot 

®v  Mass  flow  rate  per  unit  area  through 

valve 

P  Static  pressure  in  central  duct 

Pcrl  Critical  pressure  at  valve  exit 
Paxt  External  pressure  at  slot  exit 

P0  Total  pressure  in  duct 

T  Duct  static  temperature 

T«  Duct  total  temperature 

t  Time 
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Vector  of  dependent  variables 
Vector  of  diffusion  terms 
Average  duct  velocity  in  spanwlae  or  x 
direction 

Approximate  average  duct  velocity  in 

chordwlse  or  y  direction 

Slot  width 

Spanwise  coordinate 

Coordinate  at  duct  entrance 

Coordinate  at  duct  tip 

Chordwlse  coordinate 

Ratio  of  specific  heats 

Density 

Angular  frequency 


Subscripts 


Variable  from  isentroplc  flow  at  valve 
exit 

Variable  in  pressure  supply  plenum  or 
rotor  hub 


Introduction 

Circulation  control  aerodynamics  applied  to 
rotary-wing  aircraft  offers  the  advantages  of  a 
simple  and  effective  solution  to  such  problems  as 
high  vibration  levels,  retreating  blade  stall,  hub/ 
pylon  drag  and  Incorporating  higher  harmonic  con¬ 
trol.  When  the  technology  is  applied  to  stopped- 
rotor,  VTOL  aircraft  such  as  the  X-wlng,  the  his¬ 
torical  stopped-rotor  limitations  of  aeroelastlc 
divergence,  flutter,  blade  dynamic  instability,  and 
critical  resonance  conditions  during  rpm  reduction 
are  alleviated.  Reader,  et  al1  trace  the  develop¬ 
ment  of  circulation  control  rotor  technology  from 
its  inception^  through  its  application  in  the  re¬ 
cent  past  to  programs  at  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center  (DTNSRDC). 
Current  studies  at  DTNSRDC  Involve  application  of 
circulation  control  to  rotary-wing  aircraft  for  a 
variety  of  missions. 

The  concept  of  circulation  control  (CC)  as 
applied  to  helicopters  and  rotary-wing  aircraft 
which  convert  to  fixed-wing  aircraft  by  stopping 
the  rotor  are  essentially  the  aame.  Both  types  of 
aircraft  utilize  a  shaft-driven  rotor  with  blades 
having  CC  airfoils  which  generate  lift  through  the 
Coanda  effect.  The  CC  airfoils  employ  a  rounded 
trailing  edge  with  a  thin  jet  of  air  tangetlally 
ejected  from  a  spanwise  slot  adjacent  to  the  rounded 
(Coanda)  surface.  The  jet  of  air  suppresses  boundary 
layer  separation  and  moves  the  rear  stagnation 
streamline  toward  the  lower  surface,  thereby  in¬ 
creasing  lift.  Lift  is  Increased  in  proportion  to 
the  mass  flow  rate  of  compressed  air  in  the  jet. 

Pitch  and  roll  control  requirements  are  obtained 
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and  the  flux  times  area  vector  la  given  by 


by  cyclic  modulation  of  the  mess  flow  rate  with 
valves.  Higher  harmonic  cyclic  control  can  be 
similarly  applied  for  reducing  blade  stresses, 
transmitted  shears,  vibration  and  power  require¬ 
ments. 

In  view  of  the  importance  of  the  system  for 
pneumatic  control  and  distribution  of  Coanda 
blowing  air  in  a  CC  airfoil,  suprisingly  little 
attention  has  been  given  In  the  open  literature 
to  the  physics  of  these  systems.  Therefore,  to 
create  a  knowledge  base  sufficient  for  preliminary 
design  requires  basic  experimentation  aimed  at 
understanding  the  phenomenology  involved  in  the 
control  and  distribution  of  airflow  In  CC  rotors, 
l.e.,  the  "pneumodynamics"  of  CC  rotors.  The  re¬ 
sults  of  these  experiments  can  also  serve  as  base¬ 
line  data  against  which  theoretical  modals  develop¬ 
ed  for  design  analysis  can  be  evaluated.  Reader 3 
conducted  such  experiments,  in  which  mass  flow 
and  pressure  measurements  were  obtained  for  a 
simple  mock-up  of  a  rotor-blade  and  Its  pneumatic 
control  valves. 
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The  vector  of  nonhomogeneoua  tents,  which  includes 
the  effects  of  blade  rotation  at  angular  velocity 
u,  wall  skin  friction  factor  Cf(x),  wall  heat  trans¬ 
fer  coefficient  h(x),  and  mass  flux  through  the 
slot  mg(x,t)  Is 
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Theoretical  models  for  computer  simulation 
of  the  pneumodynamics  of  CC  airflow  control  sys¬ 
tems  and  rotors  have  been  developed  by  a  few 
DTNSRDC  contractors.  Perhaps  the  most  noteworthy 
among  these  are  two  computer  codes  which  include 
modeling  of  the  complex,  unsteady,  distributed 
fluid  dynamic  phenomena  In  the  rotor  blade  ducting 
system,  the  SUPERFLOWS  code  developed  by  the  Lock¬ 
heed  California  Company4  and  the  more  recent  High 
Frequency  Pneumodynamic  Analysis  (HFPA)  code  de¬ 
scribed  herein.  The  HFPA  code  represents  an  at¬ 
tempt  to  improve  on  the  physical  and  numerical 
modeling  used  in  the  earlier  SUPERFLOWS  code,  by 
eliminating  modeling  Inconsistencies  and  by  adopt¬ 
ing  an  efficient  contemporary  computational  fluid 
dynamics  formulation  suitable  for  applications 
Involving  high  frequency  cyclic  control. 
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For  completeness  the  vector  of  axial  diffusion  terms 
la  included,  although  It  is  Insignificant  In  the 
present  problem.  It  Is 


The  equation  of  state  la 


The  present  paper  gives  some  details  of  the 
physical  and  numerical  model  used  In  the  HFPA  code 
and  compares  computed  results  obtained  by  applying 
the  code  to  Reader's^  experimental  configuration 
with  the  experimental  data.  A  description  of  the 
experiment  Is  included. 


2  v  ‘ 

p  -  (y-1)  (e  -  2^  -  p  ~2~  )  <2> 

In  the  above  relations  the  conservative  variables 
m  and  e  are  defined  as 
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Basic  Equations 

The  HFPA  code  solves  the  quasl-one-dlmension- 
al,  compressible  fluid  dynamic  equations  in  con¬ 
servation  law  form,  for  unsteady  flow  In  a  span- 
wise  blowing  air  supply  duct.  Internal  to  a  rotor 
blade  with  a  Coanda  slot.  These  equations  have 
the  representation 

U+xF+C  +  V-  0  (1) 

where  the  vector  of  dependent  variables  1* 


e 


P(cvT 
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and  v,(x,t)  and  m^fx.t)  are  the  local  (with  respect 
to  spanwlsa  x  location)  slot  width  and  blowing  mass 
flux,  respectively 

The  small  average  chordwlse  velocity  component 
inside  the  duct  can  be  crudely  approximated  as 


(3) 
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In  which  HqCx)  Is  the  local  effective  height  of  the 
duct. 


The  mass  flux  distribution  through  the  Coanda 
slot  la  assumed  to  be  sufficiently  well  represented 
by  one-dimensional  isentroplc  flow  theory  with  a 
discharge  coefficient  C  (x). 


The  above  functional  relationship  is  the  usual 
expression  for  isentroplc  flow  expanding  through 
a  nozzle  from  a  plenum  (in  this  case  the  duct  in¬ 
terior)  to  ambient  conditions  (in  this  case  the 
exterior  of  the  rotor  blade).  pQ  and  Tq  are  the 
local  total  pressure  and  temperature  respectively 
and  Pext(*»0  1®  the  external  pressure  at  the  exit 
of  the  Coanda  slot.  For  choked  slot  conditions 
the  dependence  of  Eq.  (4)  on  Pgxt  i®  eliminated. 

The  blowing  slot  for  a  circulation  control 
rotor  is  generally  designed  to  be  flexible,  l.e., 
its  width  at  a  particular  spanwlse  location  can 
depend  on  Internal  duct  pressure  at  the  slot  en¬ 
trance,  external  pressure  at  Che  slot  exit  and  the 
local  pressure  distribution  over  the  rear  surface 
of  the  airfoil.  This  type  of  slot  can  easily  be 
accommodated  within  the  framework  of  the  present 
model;  the  HFPA  code,  although  used  herein  to  ob¬ 
tain  simple  constant  slot  height  results,  contains 
provisions  for  calculations  involving  the  more 
general  case  of  a  flexible  slot. 


The  computer  code  will  also  handle  the  more 
general  case  wherein  minor  flow  losses  occur  be¬ 
tween  the  air  supply  duct  and  the  slot;  for  ex¬ 
ample,  due  to  partial  obstructions  caused  by 
structural  members  interposed  between  the  air  sup¬ 
ply  duct  and  the  slot  entrance. 
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Boundary  Conditions 


At  the  upstream  boundary  x  ■  xe>  a  plenum 
supplies  Coanda  blowing  air  Co  the  duct  through  a 
valving  system.  The  valve  opening  area  has  con¬ 
stant  and  rotor-azlmuth-depcndent  components  regu¬ 
lated  by  the  collective  and  cyclic  components  re¬ 
spectively  of  the  aircraft  control  system.  The 
mass  flow  and  duct  pressures  are  thus  determined 
by  these  valve  settings.  Figure  1  shows  the  three 
different  types  of  idealized  valve  models  consider 
ed  in  the  present  analysis,  where  Ay(t)  is  the  con¬ 
trolled  value  of  the  valve  opening  and  Ae  is  the 
fixed  effective  opening  at  the  duct  entrance. 

Like  the  slot  flow,  one-dimensional  isentroplc 
flow  theory  with  a  discharge  coefficient  is  used 
to  represent  flow  through  the  valve.  The  valve 
mass  flux  is 
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where  p^  and  T  are  the  total  pressure  and 
pi  Pi 

temperature  in  the  air  supply  plenum  in  the  rotor 
hub.  For  choked  flow  through  the  valve  there  is 
no  dependence  on  the  back  pressure  in  the  duct 
entrance  downstream  of  the  valve,  p(xe,t).  The 
critical  pressure  in  the  duct  at  the  valve  for 
choked  flow  la 


Pcrit  ’  p0. 


Pi 


(  — 
\Y+lj 


y/y-1 


(6) 


An  expansion  wave  will  occur  at  the  valve  exit  for 
choked  flow.  The  pressure  pv  Just  upstream  of 
the  valve  exit  Is 


Fig.  1  Valve  Flow  Models 


pcrlt  for  p(Vt5  «  pcrit 


(8) 


The  remaining  dependent  variables  at  the  same  loca¬ 
tion  can  be  calculated  from 
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pi 


P(x.,t)  for  "<*.'*>  51  pcrlt 


(7) 


The  flux  vector,  F,  in  Eq.  (1)  must  be  given  special 
attention  at  the  duct  entrance  to  ensure  a  consistent 


■--i- 


transition  from  the  valve  to  tha  duct.  For  valve 
types  1  and  2  in  Fig.  3 


F(*e,t)  -  ayAv 


p(xa.t)Al 


For  valve  type  3  in  Fig.  1,  Eq.  (11)  aust  be  aou. 
fled  so  that  tha  second  entry  in  the  first  coluan 
vector  on  tha  right  hand  side  is  sero. 

The  downstream  boundary  condition  at  x  -  x  , 
the  rotor  tip,  aust  take  into  account  a  possible 
blowing  slot  with  area  A  incorporated  into 
the  tip  cap.  *tlp 

Therefore,  at  the  tip 

"tip  ■  V.  ■1(Po(*t.‘).To(*t.t).»wt(xft)) 

tip 


F(xt,t) 


p(«t.t) 


e(*t,t)  +  p(*t,t) 


+  p(xt.t)A(»t) 


where  A  •  0  if  there  is  no  tip  elot. 
tip 

Numerical  Analysis 


Finite-Difference  Method 

In  applying  the  theory  of  the  previous  sec¬ 
tion  to  the  analysis  of  a  rotor  blade,  the  entire 
blade,  including  the  transition  duct  (ducting  from 
hub  plenum  to  airfoil  portion  of  rotor  blade)  is 
discretized  by  dividing  it  into  spanwlse  (radial) 
segments.  The  governing  differential  equations 
are  solved  by  finite  differences  et  the  grid 
points  located  at  segment  boundaries,  including 
the  valve  exit  and  tha  blade  tip. 

Implicit  finite-difference  methods,  in  princi¬ 
ple,  offer  the  advantage  of  no  restriction  on  in¬ 
tegration  time  step.  Since  in  the  present  prob¬ 
lem  time-dependent  solutions  of  the  governing 
equations  are  of  Interest  for  dynamic  problems  of 
system  response  to  high  frequency  control  Inputs, 
and  it  is  not  desirable  to  link  integration  time 
step  to  the  spstlal  grid,  an  implicit  method  was 
chosen.  Candidate  methods  include  those  developed 
by  Briley  and  McDonald,*  learn  and  Harming,*  and 
more  recently  by  McCormack.  ?•*  In  the  preaent  re¬ 
search,  the  method  of  Beam  and  Harming*  was 


selacted  and  applied  to  the  governing  equations 
in  the  form  given  by  Eq.  (1).  In  applying  this 
method,  the  rotor  duct  Is  discretized  such  that  it 
is  divided  into  J  elements  with  J+l  grid  points,  in¬ 
cluding  a  grid  point  at  the  entrance  to  the  duct  and 
one  at  the  tip.  Following  Beam  and  Harming,”  the 
solution  is  advanced  from  time  level  n  to  time  level 
n  to  tine  level  n+1  implementing  the  "delta"  form 
of  the  difference  equations  where  at  each  grid  point 
1 
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The  contribution  of  diffusion  terms  is  assumed  minor 
in  Importance  and,  tharefore,  sufficiently  well 
approximated  by  an  explicit  evaluation  with  no  seri¬ 
ous  impact  on  atablllty.  In  Eq.  (IS),  the  implicit 
Beam  and  Harming  three-level,  second-order-eccurate 
scheme  is  obtained  by  apeclfying  6-1,  1-1/2.  A  ful¬ 
ly  implicit  one-level  scheme  suitable  for  starting 
the  calculating  or  changing  integration  tine  step  is 
obtained  by  specifying  6-1,  (-0.  Defining  the 
Jacob Ians 


If  t  “d  bi  '  If 


and  employing  modified  central  differencing  on  a 
variable  grid  where,  for  example. 


F  -F 
1+1  1-1 


then  Eq.  (IS)  leads  to  the  following  system  of  dif¬ 
ference  equations  which  will,  when  modified  using 
the  boundary  conditions,  have  block- tridiagonal  form. 
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where  L,  M,  and  N  are  the  3x3  matrices. 
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and  the  colion  vector  la 

k(fe£r)  A 
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where  numerical  damping  aa  recoaaended  by  Beam  and 
Warming6  haa  been  added  In  the  form 
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In  the  preaent  work,  after  8ame  experimentation, 
the  damping  coefficient  0  waa  aet  equal  to  zero. 

Numerical  Boundary  Conditions 

At  the  upstream  boundary,  the  boundary  condi¬ 
tions  specified  as  Eqs.  (5),  (9),  and  (10)  are  In¬ 
sufficient  for  the  numerical  solution.  Three  con¬ 
ditions  are  needed  whereas  only  two  of  these  equa¬ 
tions  are  Independent.-  A  common  and  physically  con¬ 
sistent  procedure  for  obtaining  the  extra  condition 
needed  for  flow  at  a  subsonic  Inlet  Is  to  use  one¬ 
sided  differences  as  applied  by  White  and  Anderson1* 
or  to  use  the  method  of  characteristics.  However, 

In  the  present  calculation  tha  presence  of  the  valve 
at  the  Inlet  apparently  renders  these  procedures 
useless  since  they  could  not  be  made  to  work. 

Another  approach,  which  proved  to  be  more  fruit¬ 
ful,  la  to  derive  an  approximate  auxiliary  relation¬ 
ship  based  on  Ignoring  the  time-dependent  term  In 
the  momentum  equation  between  the  first  two  grid 
points.  The  consequence  of  this  la  to  produce  an 
error  in  the  calculation  of  transients  between  the 
first  two  grid  points.  The  global  effect  of  this 
approximation  can  be  minimized  by  maintaining  a  re¬ 
latively  fine  grid  spacing  between  these  points, 
where  possible.  In  deriving  this  approximation, 
the  steady  flow  momentum  equations  corresponding  to 
the  three  valve  types  In  Fig.  1  are  written  between 
the  grid  point  1*1  at  the  entrance  and  the  next 
grid  point  1*2  Immediately  downstream.  The  result¬ 
ing  expressions  are: 


Type  1  valve 
_  2 


(p  +p2)  m 

57  CvAv  +  P1*V  +  ~2~  <VV  ’  if  A2 

(Aj+Aj)  p2u2  2  2 

*  P2A2  +  (x12-*22) 


'2  P2Dg  <VX2> 


Type  2  valve 


(P,+P2) 

■T-  Cv*v  +  P1A.  +  ~2  <W  ’  f~  ^2A2 


(A1+A2)  p/  2 

+  2  2  (xl  ~X2  )-2  p  D„(  l~x2)  (18) 


Type  3  valve 


plAe  + 


(pl'fp2) 

2 


VV  *  *2  <r  p2*2 


2  2 
PjW  2  2  Cf®2 

+  -  e-2—  (x2-x2)-  2  ^  (x1-x2) 

(19) 

In  the  above  expressions.  In  general,  •  x ,  and 
P1  "  -  pv 

where  the  inequality  holds  for  choked  flow. 

To  apply  the  upstream  boundary  conditions,  the 
differential  form  of  Eq.  (16),  (17)  or  (18)  is  com¬ 
bined  with  Eq.  (2)  and  with  the  differential  forms 
of  Eqs.  (5)  and  (10).  This  procedure  transforms 
the  boundary  conditions  into  delta  form  and  enables 
the  dependent  variables  at  1*1,  AU^  to  be  eliminated 
from  the  delta  difference  equation  for  1*2. 


i““x  tlP  Btlp 


eimax~I  *  V 


where  mtlp  Is  evaluated  from  the  total  quantities 
at  i-lmax"1,  Iad  Plaax  lx  calculated  from  an  approxi¬ 
mate  difference  form  of  the  momentum  equation 
written  between  l^g-l  and  l„r,.  The  Jacobian  term 
at  1  ,  alMX  AUi  is  taken  equal  to  its  value 

at  lmny-1  with  a  correction  for  the  difference  In 
mass  flow  between  the  two  points.  The  net  effect 
of  these  approximations  la  to  eliminate  the  de¬ 
pendent  variables  at  l"!***’  ^l«i»  from  the  delta 
difference  equation  for  l-x^x-l. 


By  the  manipulations  Just  described,  the  delta 
difference  equation  formulations  for  the  new  bounda¬ 
ry  points  1*2  and  l*lw..-l  contain  first-order  ac¬ 
curate  approximations  tor  the  physical  boundary 
conditions.  When  these  two  equations  are  combined 
with  Eq.  (15)  for  the  Interior  grid  points,  they 
constitute  an  block  trldlagonal  sat  of  difference 
equations.  A  relatively  efficient  and  simple 
direct  solution  algorithm  can  be  used  to  solve  such 
a  system  of  equations. 


W-v.W 


Initial  Conditions  and  Convergence 

The  basic  approach  followed  In  obtaining  solu¬ 
tions  to  control  inputs  consisting  of  constant 
collective  settings  or  constant  collective  settings 
with  superimposed  constant  amplitude  cyclic  inputs 
is  to  obtain  the  steady  state  or  time  periodic  solu¬ 
tions  by  allowing  them  to  evolve  from  a  transient. 

The  code  starts  with  the  blade  essentially  at 
pressure  supply  conditions.  The  initial  valve 
position  is  assumed  to  be  fully  open  as  specified 
by  its  Fourier  components  supplied  as  input.  The 
solution  is  then  advanced  with  time  as  the  flow  in 
the  blade  adjusts  itself  to  reach  steady  state  flow 
conditions  for  the  fully  open  valve.  If  cyclic 
components  are  to  be  considered,  the  cyclic  is  im¬ 
posed  after  full-open  steady  state  is  achieved. 

The  solution  is  thereafter  advanced  further  with 
time  until  periodic  conditions  are  attained.  Experi¬ 
ence  indicates  that  two  cycles  of  the  fundamental 
frequency  are  usually  sufficient  to  achieve  the 
desired  periodicity;  hence,  the  final  solution  is 
taken  to  be  the  results  during  the  third  cycle.  A 
harmonic  analysis  of  pressures,  temperacures ,  and 
mass  flow  races  is  performed  on  the  results  from 
this  third  cycle. 


The  analytical  formulation  and  the  numerical 
methods  described  in  the  previous  sections,  which 
are  employed  in  the  HFPA  code,  were  validated  by 
applying  the  code  to  a  basic  configuration  for  which 
e as lly-interpre table  experimental  results  were 
available.3  In  this  experiment  an  investigation 
was  performed  using  a  mock-up  of  a  rotor  blade  pneu¬ 
modynamic  system  incorporating  a  valve  similar  to 
the  type  3  valve  of  Fig.  1.  The  model  consisted  of 
a  supply  plenum  and  two  rotating  cams  mounted  on  the 
same  shaft  (allowing  a  mix  of  one  or  two-per-rev 
profiles  to  be  selected)  modulating  flow  from  the 
plenum  through  a  single  nozzle  into  a  stationary 
plugged  pipe  having  an  adjustable  slot.  The  pipe 
incorporated  three  total  and  two  static  pressure 
taps  distributed  along  the  length  of  the  pipe.  The 
general  arrangement  and  principal  dimensions  of  the 
model  are  shown  in  Fig.  2. 


Dynamic  measurements  of  pipe  pressure  were 
obtained  by  positioning  one  of  the  two  cams  (or 
portions  of  both)  over  the  nozzle  inlet,  setting  a 
plenum  pressure  and  obtaining  data  for  a  range  of 
cam  rotational  speeds.  In  addition,  quasi-dynamic 
results  were  obtained  by  turning  the  cam  in  discrete 
30-deg  increments,  each  tlaa  recording  data  for 
this  non-rotating  condition.  Average  mass  flow  data 
were  collected  from  a  venturi  located  upstream  of 
the  plenum. 

Fig.  3  shows  the  effective  valve  area  formed 
by  the  one-  and  two-per-rev  cams  spaced  at  a  distance 
of  0.01  in.  from  the  nozzle  inlet.  The  points  shown 
were  calculated  from  the  measured  cam  profile;  the 
area  variations  are  intended  to  approximate  sinusoids 
as  shown. 


«• 


flow  losses  which  occur  over  a  range  of  mass  flow 
rates  In  a  dynamic  or  quasl-dynamlc  analysis. 


after  in  elapsed  tine  of  0.025  seconds  and  the 
tlmesiup  was  reduced  to  1200  steps  per  revolution 
thereafter.  Execution  consumed  approximately  two 
and  one-half  minutes  IBM  3033S  CPU  time.  In  spite 
of  the  implicit  finite-difference  scheme,  larger 
tlmesteps  led  to  numerical  Instability.  Prelimi¬ 
nary  calculations  made  to  compare  HPPA  and  SUPER- 
FLOWS  results  for  the  sample  case  of  Ref.  4  indi¬ 
cate  that,  although  the  present  implicit  method  is 
still  restricted  to  small  tlmesteps,  these  time- 
steps  are  several  times  larger  than  those  permitted 
by  the  multistep  explicit  method  used  in  SUPER¬ 
FLOWS. 

Typical  quasl-dynamlc  numerical  and  experi¬ 
mental  results  indicating  the  dependency  of  mass 
flow  rates  on  nozzle  control  area  at  a  constant 
slot  height  of  0.042  in.  are  displayed  for  the 
one-per-rev  cam  in  Fig.  4  and  the  two-per-rev  cam 
in  Fig.  S.  The  mass  flow  rates  are  plotted  against 
the  discrete  azimuthal  cam  settings  at  which  the 
readings  and  calculations  were  made.  As  indicated 
on  the  figures  the  plenum  pressures  were  set  to 
different  values  for  the  two  cams.  The  figures 
indicate  the  ability  of  the  cam  control  valve  to 
modulate  mass  flow  so  as  to  track  the  cyclic  vari¬ 
ation  in  valve  setting. 
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4  Quasi-Dynamic  Mass  Flow  Rate  vs  Cam 
Position  for  Ona-Per-Rev  Cam 


Figs.  4  and  5  also  show  tha  effect  of  as- 
summed  constant  values  of  0.95  and  0.75  for  equal 
slot  and  valve  discharge  coefficients  in  the  KFPA 
code.  The  results  obtained  for  the  higher  values 
Cv  ■  C§  •  0.95  agree  more  favorably  with  the  ex¬ 
perimental  results,  despite  the  fact  that  the 
lower  discharge  coefficient  is  likely  more  physi¬ 
cally  realistic.  Strict  Interpretation  of  C,  and 
Cy  as  discharge  coefficients  is  probably  not  ds- 
sirabla  since  the  various  deficiencies  and  approxi¬ 
mations  in  the  computer  model  can,  to  some  degree, 
be  absorbed  by  adjusting  them.  Moreover,  the  con¬ 
cept  of  a  constant  nostle  discharge  coefficient  is 
at  best  only  a  crude  approximation  for  the  actual 
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Fig.  5  Quasi-Dynamic  Mass  Flow  Rate  vs  Cam 
Position  for  Two-Per-Rev  Cam 

Figs.  6  and  7  show  the  quasl-dynamlc  results 
for  total  pressure  corresponding  to  the  mass  flow 
rates  of  Figs.  5  and  6  again  plotted  against  the 
discrete  azimuthal  cam  setting  at  which  the  readings 
and  calculations  were  made.  These  plotB  also  show 
the  better  agreement  obtained  for  Cs  »  0.95  and  the 
ability  of  the  cam  arrangement  to  modulate  total 
pressure  so  as  to  track  the  cyclic  variation  in 
valve  setting.  The  total  pressures  shown  in  the 
figure  are  for  the  furthest  outboard  pressure  tap 
location;  the  measured  and  computed  variation  in 
total  pressure  along  the  length  of  the  pipe  was 
minimal. 
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Fig.  6  Quasi-Dynamic  Total  Pressure  vs  Cam 
Position  for  One-Per-Rev  Cam 

Fig.  8  is  the  dynamic  counterpart  of  the 
quasi-steady  total  pressures  plotted  in  Fig.  6 
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and  Is  typical  of  the  dynamic  results  obtained. 
Extensive  results  of  the  dynamic  calculations  and 
measurements  will  be  reported  In  a  subsequent 
paper;  space  limitations  will  not  permit  their 
presentation  here.  The  numerical  results  of  Fig. 
8  were  obtained  with  Cy  «  C«  ■  0.95.  The  numeri¬ 
cally  computed  pressure  profiles  are  actually 
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Quasi-Dynamic  Total  Pressure  vs  Cam 
Position  for  Two-Per-Rev  Cam 


displaced  laterally  in  azimuth  position  from  the 
valve  area  profiles  of  Fig.  3  due  to  an  approxi¬ 
mately  35  deg  computed  phase  shift.  This  phase 
shift  is  principally  due  to  the  "capacitance  lag" 
effect  caused  by  the  finite  pipe  volume.  A 
smaller  portion  of  the  delay  is  due  to  the  finite 
speed  of  wave  propogatlon,  so-called  "sonic  lag." 
The  phase  difference  between  the  valve  area  set¬ 
ting  and  the  pressure  response  to  it  has  been  eli¬ 
minated  from  Fig.  8  for  purposes  of  comparison 
with  the  experimental  data.  Although  present  in 
the  experiment,  phase  shift  was  not  an  objective 
of  the  experimental  measurements  and  is  not  availa¬ 
ble  from  the  data.  In  Figs.  9  and  10  the  quasi- 
dynamic  results  of  Fig.  6  are  superimposed  onto 
the  dynamic  results  to  demonstrate  the  effect  of 
unsteady  flow  on  the  magnitude  and  shape  of  the 
experimental  and  numerical  pressure  responses, 
respectively.  The  unsteady  effect  seems  to  cause 
a  slight  decrease  in  the  peak  pressure  from  the 
quasi-dynamic  results  in  both  the  numerical  and 
experimental  data.  However,  this  decrease  is 
fairly  small;  hence,  except  for  the  phase  dif¬ 
ference,  the  quasi-dynamic  results  approximate  the 
unsteady  results  reasonably  well. 

Conclusions 
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Fig.  8  Dynamic  Cyclic  Variation  of  Total 
Pressure 
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Comparison  of  Experimental  Dynamic  and 
Quasi-Dynamic  Pressure  Profiles 


The  results  presented  herein  tend  to  indicate 
that  quasl-one-dimensional  flow  theory  can  be 
applied  to  predict  with  reasonable  accuracy  the 
pneumodynamics  of  a  circulation  control  rotor. 

The  implicit  finite  difference  technique  used  re¬ 
presents  an  improvement  over  an  explicit  method 
used  in  a  previous  computer  model  but  is  still 
restricted  to  small  timesteps.  The  results  further 
indicate  that  simple  cam- type  valving  can  be  used 
to  modulate  the  flow  and  pressure  of  Coanda  blow¬ 
ing  air  to  approximately  track  the  variation  in 
the  valve  opening  area.  Quasi-dynamic  numerical 
and  experimental  results  can  predict  dynamic  re¬ 
sults  to  a  good  approximation,  but  do  not  reflect 
the  effects  of  capacitance  and  sonic  lag. 
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Fig.  10  Comparison  of  Numerical  Dynamic  and 
Quasi-Dynamic  Pressure  Profiles 
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